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Abstract 

We present a new approach to study the properties of the sun. We consider small variations of the physical 
and chemical properties of the sun with respect to Standard Solar Model predictions and we linearize the 
structure equations to relate them to the properties of the solar plasma. By assuming that the (variation 
of) the present solar composition can be estimated from the (variation of) the nuclear reaction rates and 
elemental diffusion efficiency in the present sun, we obtain a linear system of ordinary differential equations 
which can be used to calculate the response of the sun to an arbitrary modification of the input parameters 
(opacity, cross sections, etc.). This new approach is intended to be a complement to the traditional methods 
for solar model calculation and allows to investigate in a more efficient and transparent way the role of 
parameters and assumptions in solar model construction. We verify that these Linear Solar Models recover 
the predictions of the traditional solar models with an high level of accuracy. 



1 Introduction 



In the last three decades, there was an enormous progress in our understanding of the sun. The 
predictions of the Standard Solar Model (SSM), which is the fundamental theoretical tool to investigate 
the solar interior, have been tested by solar neutrino experiments and by helioseismology. 

The deficit of the observed solar neutrino fluxes, reported initially by Homestake and then 
confirmed by GALLEX j2] and SAGE |3| (which subsequently merged into GNO jl]), Kamiokande 
[5] and Super-Kamiokande [6], generated the so-called "solar neutrino problem" which stimulated a 
deep investigation of the solar structure (see e.g. [7]). The problem was solved in 2002 when the 
SNO experiment [8] obtained a direct evidence for flavour oscillations of solar neutrinos and, moreover, 
confirmed the SSM prediction of the neutrino flux with an accuracy which, according to the latest 
data [9], is equal to about ~ 6% (see e.g. p!0]l[^ 

At the same time, helioseismic observations have allowed to determine precisely several important 
properties of the sun, such as the depth of the convective envelope which is known at the ~ 0.2% level, 
the surface helium abundance which is obtained at the ~ 1.5% level and the sound speed profile which is 
determined with an accuracy equal to ~ 0.1% in a large part of the sun (see e.g. [161 EZ] and references 
therein). As a results of these observations, the solar structure is now very well constrained, so that 
the sun can be used as a solid benchmark for stellar evolution and as a "laboratory" for fundamental 
physics (see e.g. [B]). 

The future times could be even more interesting. KamLAND reactor (anti)neutrino experiment 
|19] has confirmed the flavour oscillations hypothesis and has refined the determination of neutrino 
parameters. We now reliably know the solar neutrino oscillation probability and we can go back to 
the original program of solar neutrino studies, i.e. to probe nuclear reactions in the solar core. The 
present and future solar neutrino experiments, such as Borexino |2^ and SNO-I- |21], have the potential 
to provide the first direct measurements of the CNO and pep neutrinos, thus probing dominant and 
sub-dominant energy generation mechanisms in the sun. 

At the same time, a new solar problem has emerged. Recent determinations of the photospheric 
heavy element abundances |221 123] indicate that the sun metallicity is lower than previously assumed 
|24] . Solar models that incorporate these lower abundances are no more able to reproduce the helio- 
seismic results. As an example, the sound speed predicted by SSMs at the bottom of the convective 
envelope disagrees at the ~ 1% level with the value inferred by helioseismic data (see e.g.j25j). Detailed 
studies have been done to resolve this controversy (see e.g.[T71 |23]), but a definitive solution of the 
"solar composition problem" still has to be obtained. 

In this framework, it is important to analyze of the role of physical inputs and of the standard 
assumptions for SSM calculations. This task is not always possible in simple and clear terms, since 
SSM construction relies on (time-consuming) numerical integration of a non linear system of partial 
differential equations. Several input parameters are necessary to fully describe the property of the solar 
plasma and some of them are not single numbers but complicated functions (like e.g. the opacity of the 
solar interior) which, in principle, can be modified in a non trivial way. Moreover, any modification of 
the sun produces a variety of correlated effects which have to be taken into account all together if we 
want to correctly extract information from the comparison of theoretical predictions with observational 

^ For the sake of precision, the first model-independent evidence for solar neutrino oscillations and the first determination 
of the *B solar neutrino flux has been obtained in 2001 (see e.g. [TT]) by comparing the SNO charged-current result [T2] 
with the SK data with the method proposed by [T3]. The year 2002 is, however, recognized as the 'annus mirabilis' [TJ] 
for the solar neutrino physics. During 2002, in fact, the SNO neutral-current measurement El and the first KamLAND 
results [12] were released and, moreover, Raymond Davis, Jr. and Masatoshi Koshiba were awarded with the Nobel Prize. 
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data. 

In order to overcome these difficulties, we provide a tool which is at same time simple and accurate 
enough to describe the effects of a generic (small) modification of the physical inputs. The starting 
point is the fact that, despite the present disagreement with helioseismic data, the SSM is a rather 
good approximation of the real sun. We can thus assume small variations of the physical and chemical 
properties of the sun with respect to the SSM predictions and use a linear theory to relate them to the 
properties of the solar plasma. With the additional assumption that the (variation of) the present solar 
composition can be estimated from the (variation of) the nuclear reaction rate and elemental diffusion 
efficiency in the present sun, we obtain a linear system of ordinary differential equations, which can be 
easily integrated. 

We believe that the proposed approach can be useful in several respects. First, the construction of 
these Linear Solar Models (LSM) can complement the traditional methods for solar model calculations, 
allowing to investigate in a more efficient and transparent way the role of parameters and assumptions. 
In a separated paper, we will use LSMs to discuss in general terms the role of opacity and metals in 
the solar interior [26]. Moreover, it can help to introduce new effects and to understand their relevance, 
prior to the implementation into the more complicated SSM machinery. Finally, this simplified approach 
could open the field of solar physics beyond the small community of the SSM's builders with the result of 
making the sun a more accessible "laboratory" . Clearly, LSMs have not to be intended as an alternative 
to SSMs, which remain the fundamental theoretical tool to compare with observations. 

The plan of the paper is the following. In the next section we briefly review SSM calculations. In 
sect J3] we expand to linear order the structure equations of the present sun. In sect j4] we calculate the 
properties of the sun and of the solar plasma that are necessary to define LSMs. In sect jS] we give the 
integration conditions for the linearized equations. In sectjH] we discuss how to estimate the chemical 
composition of the sun. In sect J7] we give the equations that define LSMs in their final form. Finally, 
in sect js] we compare the results of LSMs with those obtained by using the standard method for solar 
model calculations, showing the validity of our simplified approach. 



2 The standard solar model 

In the assumption of spherical symmetry, the behavior of the pressure (P), density (/)), temperature 
(T), luminosity (/) and mass (m) in the sun is described by the structure equations |27] : 



dm 
dr 

dP _ _G^m 



-p 



m 

dr 

dT G^mTp 



47rrV e(p,r,Xi 



V 



dr r^P 
P = P{p,T,Xi) (1) 

where r indicates the distance from the center, e{p,T,Xi) is the energy produced per unit time and 
mas^ the function P = P{p, T, Xi) describes the equation of state (EOS) of solar matter and Xi{r) are 

^We neglect, here and in the following, the contribution to energy balance due to the heat released by the various shells 
of the sun in their thermodynamical transformations. This contribution, which depends on the time derivative of pressure 
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the mass abundances of the various chemical elements inside the sun. In the equation which describes 
the energy transport, the temperature gradient V is defined as V = dlnT/dln P. If the energy transfer 
is due to radiative processes (i.e. at the center of the sun) one has: 



V = Vrad 



(2) 



where the radiative gradient Vrad is given by: 



Vrad 



IGvracGN 
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K{p,T,Xi)lP 



(3) 



and n{p,T, Xi) is the opacity of solar matter. In the presence of convective motions (i.e. in the outer 
layer of the sun), the value of V has to be calculated by taking into account the convective energy 
transport which generally provides the dominant contribution. In a large part of the sun convective 
envelope, one has 



This is expected when convection is very efficient and a negligible excess of V over the adiabatic value 
Vad is sufficient to transport the whole luminosity. In the outermost layers of the sun, the situation 
is more complicated. The precise description of convection in this regime is still an unsolved problem. 
One uses a phenomenological model which predicts the efficiency of convection as a function of the 
"mixing length parameter" a, which is related to the distances over which a moving unit of gas can be 
identified before it mixes appreciably (see e.g. I2'7j). The transition between the internal radiative core 
and the external convective envelope occurs at the radius i?b where: 



as it is prescribed by the Schwarzschild criterion which states that convective motions occur in the 
region of the star where Vrad('') > Vad- 

The chemical composition of the present sun is not known but has to be calculated by coupling 
eqs. ([T|, which describe the mechanical and thermal structure, with the equations that describe the 
chemical evolution of the sun (see e.g. [27j for details). It is generally assumed that the sun was 
born chemically homogeneous and has modified its chemical composition due to nuclear reactions and 
elemental diffusion. One assumes that relative heavy element abundances in the photosphere have not 
been modified during the evolution and uses the observationally determined photospheric composition 
(see e.g.|22| |23]) to fix the initial heavy element admixture. We remark that the chemical abundances 
are uniform in the convective region due to the very efficient mixing induced by convective motions. 
The values Xi^^, evaluated at the bottom of the convective region, are, thus, representative for the solar 
surface composition. 

If a complete information about the initial composition was available and if a complete theory 
of convection was known, then there would be no free parameter. In practice, one has three free 
parameters, namely the initial metal abundance Zjni, the initial helium abundance l^ni and the mixing 
length parameter a. These are tuned in order to reproduce, at the solar age = 4.57 Gyr [29 , the 
observed solar luminosity Lq = 3.8418 • lO^^erg/s [30], the observed solar radius Rq = 6.9598 • lO^'^cm 
|31| and the surface metal-to-hydrogen ratio (Z/X)i,. In this paper we use the AS05 composition |22] 
which corresponds to (Z/X)b = 0.0165. 

and temperature, is sub-dominant as can be understood from the fact that the sun evolves on time scales much larger 
than the Kelvin-Helmotz time. 



V ~ 



Vad - 0.4 . 



(4) 




(5) 
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A standard solar model (SSM), according to the definition of [28], is a solution of the above problem 
which reproduces, within uncertainties, the observed properties of the sun, by adopting physical and 
chemical inputs chosen within their range of uncertainties. In this paper, we refer to the SSM obtained 
by using FRANEC code [32], including: Livermore 2006 equation of state (EOS) [33]; Livermore 
radiative opacity tables (OPAL) [3l] calculated for the AS05 chemical composition; molecular opacities 
from ref. [35] and conductive opacities from ref. [36] both calculated for AS05 composition; nuclear 
reaction rates from NACRE compilation [38] , taking into account the recent revision of the astrophysical 
factors 5i^i4, 5*3^4 and S'i,7[39l HOI HI] and the ^Be electron capture rate from [42] . 



3 Linear expansion of the structure equations 

A variation of the input parameters (EOS, opacities, cross sections, etc.) with respect to the standard 
assumptions produces a solar model which deviates from SSM predictions. The physical and chemical 
properties of the "perturbed" sun can be described according to: 

h{r) = h{r)[l + 6h{r)] 
Xi{r) = Xi{r)[l + 6X,{r)] 

Y{r) = F(r) + AY{r) (6) 

where h = I, m, T, P, p and we use, here and in the following, the notation Q to indicate the SSM 
prediction for the generic quantity Q. For convenience, the deviations from SSM are expressed in terms 
of relative variations 6h{r) and 6Xi(r) for all quantities except the helium abundance, for which it is 
more natural to use the absolute variation AY(r). 

In this paper, we will be mainly concerned with modifications of the radiative opacity, n(p,T, Xi), 
and of the energy generation coefficient, e{p,T, Xi). We indicate with 5K,(r) and 6e{r) the relative 
variations of these quantities along the SSM profile, i.e.: 

I{T{r),p{r),Y{r),Xi{r))=I{T{r),p{r),Y{r),X,{r^ [1 + 61{r)] (7) 

where I = k, e. When we consider the effect of a perturbation 61{r) on the sun, we have to take into 
account that the perturbed solar model, due to that modification, has a different density, temperature 
and chemical composition with respect to the SSM. The total difference (57*°* (r) between the perturbed 
sun and the SSM at a given point r is defined by: 

IiT{r),pir),Yir),Xi{r))=IiTir),p{r),Y{r),Xi{r))[l + 5l'"\r)]. (8) 

If we consider small perturbations {61 <C /), we can expand to first order in 6T{r), Sp{r), 5Xi(r) and 
Ay(r), obtaining: 

5I^°\r) = Irir) 5T{r) + Ip{r) 5p{r) + Iy(r) Ay(r) + ^ /^(r) 5Xi{r) + 51{r) (9) 

i 

Here, the quantities Ij describe the dependence of the properties of the stellar plasma from temperature, 
density and chemical composition and are given by: 

^^""^ ainT SSM 
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dY 



SSM 



(10) 



where I = e, k, P. The symbol |ssM indicates that we calculate the derivatives Ij along the density, 
temperature and chemical composition profiles predicted by the SSM. 

By using the above notations, we linearize the structure eqs. ([T]), obtaining: 

d6m 



dr 
d6P 
dr 
5P 

d5l 
dr 

d6T 

dr 
d6T 
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[5p — 5m] 
[6m + 5p- 5P] 



1 

1 

[Pp 5p + Pt6T + PyAY + Y.Pi 5Xi\ 

i 

1 
k 
1 

h 



[(1 + ep)5p + eT5T + eyAY + ^ 6Xi -51 + (5e] 

i 

[51 + {kt - 4)(5r + {Kp + l)5p + kyAY + 5Xi + 5k] 



— [5m + 5p- 5P] 
It 



Rad. 
Conv. 



fill 



-1 



where Ih = dln{h) / dr^ represents the scale height of the physical parameter h in the SSM. The last 
two equations correspond to the linear expansion of the energy transport equation in the radiative and 
in the convective regime respectivel}|^ In our approach, we use the radiative transport equation for 
r < i2b and the convective transport equation for r > R\,, where i?b = 0.730 -R© is the lower radius 
of the convective envelope predicted by our SSM. We calculate a posteriori the relative variation 5R\j 
of the extension of the convective region by applying the Schwarzchild criterion to the solutions of 



eqs. (11). By expanding to first order, we obtain: 



5Rb 



5Vrad,h _ 5Ktot,h + 5Ph + 5h - 4(5rb - 5mb 



Cb 



Cb 



(12) 



where Cb = din Vrad(^b)/c?ln r = 11.71. Here and in the following, the subscript "b" indicates that 
the various quantities are evaluated at r = i2b) i-e. at the bottom of the SSM's convective envelope. 



4 Properties of the sun and of the solar matter 

In order to define the linearized structure equations, we have to calculate the functions Z/i(r) and the 
logarithmic derivatives Ij{r). In Figjl] we show the inverse scale heights l/lh{r) as a function of the 
solar radius. The plotted results have been obtained numerically and refer to our SSM. We see that the 
functions l//p(r) and l/lT^r) vanish at the center of the sun while they grow considerably (in modulus) 
in the outer regions, as a result of the fast decrease of pressure and temperature close to surface. The 
kink in the temperature scale height at Ry, = O.73Oi?0 marks the transition from the internal radiative 

■^In the presence of convection, we assume that V = Vad and that the adiabatic gradient is not affected by the performed 
modifications of the input parameters. 
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Figure 1: The inverse scale heights l/lh{r) of temperature (black), pressure (red), mass (green) and luminosity 
(blue) as a Junction of the solar radius. 



region to the outer convective envelope. The functions l/lm{r) and l/li{r) have the opposite behaviour; 
they diverge at the solar center and vanish at large radii, as it is expected by considering that most of 
the mass and of the energy generation in the sun is concentrated close to the solar center. To be more 
quantitative, the solar luminosity is produced in the inner radiative core (r < 0.3 -Rq), which contains 
approximately 60% of the total solar mass. The radiative region (r < i?b) which covers about 40% of 
the total volume of the sun, includes approximately 98% of the total solar mass. 

In Figj2] we show the logarithmic derivatives Ij{r) calculated numerically along the temperature, 
density and chemical composition profiles predicted by our SSM. The left panel refers to the opacity 
derivatives K,j{r) and shows that opacity is a decreasing function of temperature and helium abundance, 
while it is an increasing function of density. The coefficient kz{t) = din K/dln Z\ssm quantifies the 
dependence of opacity on the total metallicity Z. It has been calculated by rescaling all the heavy 
element abundances by a constant factor, so that the metal admixture remains fixed. Metals provide 
about ~ 40% of the opacity at the center of the sun, while they account for about 80% of the total 
opacity at the bottom of the convective region (see also [17] ) . 

In the right panel of Figj2j we present the logarithmic derivatives ej(r) of the energy generation 
coefficient. These have been calculated by assuming that the abundances of the secondary elements 
for PP-chain and CN-cycle, namely ^He, ^^C and ^^N can be estimated as it is described in Appendix 
A. The bumps at r ~ 0.28 i?0 are related to out-of-equilibrium behaviour of the ^He abundance. The 



presence of these features has, however, a negligible influence on the solutions of eqs.(ll ). The function 
l///(r), which multiplies all the ej{r) coefficients, drops, in fact, rapidly to zero at r ~ O.25i?0 

In the more internal region (where ^He assumes the equilibrium value) , the displayed results can be 



For most applications, one could assume that secondary elements have the equilibrium abundances all over the sun. 
The errors implied by this assumption are usually negligible. In the region where it is not valid, in fact, the nuclear 



reaction rates are small and, thus, a negligible amount of energy and neutrinos are produced. In some cases (see sect 8.21, 
however, it is necessary to improve this approximation. For this reason, we developed the method described in Appendix 
A, which allows to describe with good accuracy the behaviour of secondary elements all over the sun. 
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Figure 2: Left panel: The logarithmic derivative Kj{r) of the radiative opacity with respect to temperature (black), 
density (red), helium (green) and metals (blue) as a function of the solar radius. Right panel: The logarithmic 
derivative ej{r) of the energy generation coefficient with respect to temperature (black), density (red), helium 
(green) and metals (blue - multiplied by a factor 100) as a function of the solar radius. 



understood by considering that energy is produced by the PP-chain (~ 99%) with a small contribution 
(~ 1%) by the CN-cycle. The slight increase of exir) at r < 0.15 Rq reflects the fact that the CN-cycle, 
which strongly depends on temperature, gives a non-negligible contribution only at the center of the 
sun. The coefficient ep{r) is approximately equal to one, as it is expected by considering that the rate 
of two-body reactions in the unit mass is proportional to the density. Finally, the behavior of eY(r) and 
ez{r) is understood by considering that the pp- reaction rate is proportional to the hydrogen abundance 
squared, i.e. to = (1 - y - Z)^. 

We have also evaluated numerically the logarithmic derivative of the pressure Pj{r) by using the 
Livermore 2006 EOS [33]. The obtained results can be approximated with good accuracy (about 1%) by 
using the perfect gas EOS: P = (pk-QT) / m^j) , where is the atomic mass unit, k-Q is the Boltzmann 
constant and /i = (2 - 5/4 T - 3/2 Z)" Ms the mean molecular weight. One obtains: 



PT{r) 



dY 



8 - 5Y{r) - 6Z{r) 



(13) 



The derivative with respect to the metal abundance Pz{r) is of the order of few per cent and will be 
neglected in the following. By using the above relations, we obtain the following equation: 



6p{r) = 5P{r) - 5T{r) - PyAYi^r) 



(14) 



that can be used to eliminate the quantity 6p from the linear structure eqs.(ll ). We arrive at: 

1 



d6m 

dr 
ddP 

dr 
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[5P - 5T -dm-Py Ay] 
[-5T + 5m-Py Ay] 
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7 r/ -| 

— = -[(3pSP + Pt ST -61 + Py Ay + PzSZ + 6e] 
dr k 

d3T 1 

— — = — [ap5P + axST + 51 + ay l^Y + az5Z + 5k] Rad. 
dr It 

^ = ^ [-6T + 5m - Py AY] Conv. (15) 

dr It 

where the factors Oh and /3/j are given by: 

ap = Kp + 1 ut = kt — Kp — 5 ay = —{^p + l)Py + Ky az = kz 

Pp = ep + 1 PT = eT-ep-l /3y = -(e^ + l)Py + ey Pz = €z (16) 

In the above equations, we indicate with 6Z the relative variation of the total metallicity, defined by: 

Z(r) = Z(r) [1 + 6Z{r)] (17) 

and we implicitly assume that the heavy element admixture is unchanged]^ 

5 Boundary conditions 

In order to solve the linear structure equations, one has to specify the integration conditions. This can 
be done quite easily at the center of the sun. From the definitions of mass and luminosity, we know 
that m(0) = and /(O) = 0. This implies that: 



h{r) 



dh{0) 



1 + 5h(0) = lim = ^f- (18) 
^ ^ - .0 h{r) dm ^ ^ 

dr 



for h = m, L so that one obtains: 



MO) = 5Po-5To-Py,oAYo 

5P{0) = 5Po 

6T{0) = 5To 

51(0) = Pp,o5Po + PTfiSTo + Py,oAYo + Pz,o5Zo + 6eo (19) 

where the subscript "0" indicates that the various quantities are evaluated at the center of the sun. 

To implement surface conditions, we exploit the fact that the sun has a convective envelope where 
the solution of the linearized structure equations can be explicitly calculated, as it is explained in the 
following. 

First, we take advantage of the fact that there are no energy producing processes and that a 
negligible fraction of the solar mass is contained in the convective region. This implies that l///(r) = 
and l/lm{r) ~ 0, so that we can integrate the energy generation and the continuity equation. We 
obtain 5l{r) = and 5m(r) ~ where we clearly considered that only solutions that reproduce the 
observed solar luminosity and solar mass are acceptable. 

Then, we consider that the chemical abundances are constant due to the very efficient convective 
mixing, so that Ay(r) = Ayb and Py(r)Alb = —SfJ'h, where (5/ib is the relative variation of the mean 



^The effects of changes in the heavy element admixture can be described in our approach by proper variations 5k{r) 
and Se{r) (see [25]). 
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molecular weight at the bottom of the convective region. By taking this into account, we can rewrite 
the transport equation in the form: 

where 5u is the fractional variation of squared isothermal sound speed u = P/p = k^T / fim^. The 
general solution of this equation is 6u{r) = duy, [rb/T(r)] which shows that, in order to avoid 6u{r) to 
"explode" at the surface of the sun, it necessarily holds 5u{r) = 6T{r) — S/i^ — in the internal layers 
of the convective envelope. 

By using this result into the hydrostatic equilibrium equation, we obtain 5P{r) = 5p{r) ~ 5C where 
5C is an arbitrary constant. Finally, by considering that 5p{r) is approximately constant in the internal 
layers of the convective region (where most of the mass of the convective envelope is contained), we 
use the continuity equation to improve the condition 5m{r) 2± 0, obtaining 5m\) = — mconv where 
rficoTw = Mcoiiv/Mq = 0.0192 is the fraction of solar mass contained in the convective region, as resulting 
from our SSM. 

In conclusion, we have: 

(5m(i?b) = -^conv SC 

SP(R^) = SC 

ST(R^) = ,5^b = -Py,bAyb 

Sl(Rb) = (21) 

The factor 6C is a free parameter that cannot be fixed from first principles, since we cannot model 
exactly convection in the outermost super-adiabatic region. It has basically the same role as mixing 
length parameter in SSM construction. 



6 The chemical composition of the sun 

The chemical composition of the perturbed sun should be calculated by integrating the perturbed struc- 
ture and chemical-evolution equations starting from an ad-hoc chemical homogeneous ZAMS model. 
However, this would lead to several complications in our simplified approach. We prefer to use a 
simple approximate procedure that allows to estimate with sufficient accuracy the helium and metal 
abundances of the modified sun, without requiring to follow explicitly its time-evolution. 

6.1 Notations 

In order to quantify the relevance of the different mechanisms determining the present composition of 
the sun, we express the helium and metal abundance according to: 

Y{r) = Yini [1 + Drir)] + Fnuc(r) 

Z{r) = Zini [1 + Dzir)] (22) 

Here, lini and Zi^i are the initial values which, as explained previously, are free parameters in solar 
model construction and are adjusted in order to reproduce the observed solar properties. In our SSM, 
we have Y[^i = 0.2611 and Zini = 0.0140. The terms Dy^r) and Dz{r) describe the effects of elemental 
diffusion. Finally, Y^-adr) represents the total amount of helium produced in the shell r by nuclear 
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Figure 3: Left panel: The black line shows the amount of helium produced by nuclear reactions Ynuc{r)- The red 
line shows the energy generation coefficient e(r) suitably rescaled to superimpose with Ynuc{r). Right panel: The 
terms Dyir) and Dz{r) which describe the effects of elemental diffusion in the SSM. 



processes and can be calculated by integrating the rates of the helium-producing reactions during the 
sun history. 

In Figjsjwe show the quantities l^nuc('')) -Dy(r) and Dz{r) estimated from our SSM. Nuclear he- 
lium production is relevant in the internal radiative core (r < O.SRq) where it is responsible for an 
enhancement of the helium abundance which can be as large as l^nuc,o — 0.35 at the center of the sun. 
Elemental diffusion accounts for a ~ 10% increase of helium and metals in the central regions with 
respect to external layers of the radiative region. The convective envelope, due to the very efficient 
convective mixing, is chemically homogeneous and can be fully described in terms of two numbers: the 
surface helium and metal abundances, indicated with Yb and Zb respectivel}]^ These are related to 
initial values by 

lb = Yini [1 + Dy,b] (23) 
Zh = Zini [1 + Dz,h] 

where Dy^h and Dz.h parametrize the effects of elemental diffusion. In our SSM, we have Y\y = 0.229, 
Zb = 0.0125, DY,h = -0.121, Dz,h = -0.105. 

We are interested in describing how the chemical composition is modified when we perturb the SSM. 
In the radiative core (r < R^), we neglect the possible variations of the diffusion term^and we write: 



AY(r) = AYr 



mi 



l + Dvir) 



6Z{r) = 5Z,^, (24) 

where AYnuc('^) is the absolute variation of the amount of helium produced by nuclear reactions. A 
better accuracy is required in the convective region, because the surface helium abundance lb is an 



^ We use the subscript "b" to emphasize that the surface abundances coincide with the values at the bottom of the 
convective region. 

^The diffusion terms Dyir) and Dz{r) are at the few per cent level in the radiative region. Their variations are, thus, 
expected to produce very small effects on the solar composition. 
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observable quantity. We, thus, discuss explicitly the role of diffusion and we write: 

5Z^ = SZi^i + -^^^6Dz,^ (25) 

where 5-Dy,b and 5Dz,h are the fractional variations of the diffusion terms. 

It is important to remark that Alb and 5Z\y are related among each other, since the metals-to- 
hydrogen ratio at the surface of the sun is observationally fixed. By imposing (5(Z/X)b = 0, we obtain: 

5Z^ = ^AFb (26) 

l-Fb 

where we considered that Xb ~ 1 — lb- This relation can be rewritten in terms of the initial helium 
and metal abundances, obtaining: 

6Zi^i = Qo Alini + Qi 5Dy,c + Q2 SDz,^ (27) 

where the coefficients Qi are given by: 

Qi = - ^'"^^■'^ =+0.041 
1^ Fb 

Q2 = - ^fg-" =+0.118. (28) 
6.2 Production of helium by nuclear reactions 

In order to predict the helium abundance in the radiative region, one has to estimate the variation of 



nuclear production of helium Al^uc('^)) see eq.(24). We note that the quantity l"nuc(?') varies propor- 
tionally to e(r) in the SSM, as it is seen from FigjSj This is natural since the helium production rate is 
directly proportional to the energy generation rate and, moreover, the energy generation profile e(r) is 
nearly constant during the past history of the sun. We assume that the observed proportionality holds 
true also in modified solar models, obtaining as a consequence the following relation: 

Ay„uc(r-) = F„uc(r)5e*°*(r) (29) 



If we expand Se^°^{r) according to rel.([9]) and solve with respect to AY{r), we can recast eq.(24) in the 
form 13 

AY{r) = ey(r) Al^„i + ^p{r) 5P{r) + er(r) dT{r) + C.(r) 6e{r) (30) 

where: 

^ ' l-l^nuc(r)[ey(r)-Py(r)ep(r)] 



*In the expansion of 5etot{r) we neglected the term ez{r)SZ{r). This term is expected to give a negligible contribution 
since the coefficient ez(r) is at most equal to -0.05, whereas the coefficients erir), €p{r) and €Y{r) are order unity or 
larger, see right panel of figj2] 
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Figure 4: The factors ^j{r) defined in eqs.(31) 



' l-l^nuc(r)[ey(r)-Py(r)ep(r)] 

c / N ^ Ynucir) [eT{r) - ep(r)] 

^ 1-F„,c(r)[6y(r)-Py(r)ep(r)] 

^ ' l-i"nuc(r)[ey(r)-Py(r)ep(r)] ' ^ ' 

The factors '^^(r) are shown in Fig|4j 

6.3 Elemental diffusion in the convective region 

The total mass Mconv.i of the i— th element contained in convective envelope evolves due to elemental 
diffusion according to relation 



1 dM, 



conv,i 

~ H 



(32) 



Here, the parameter H is the "effective thickness" of the convective region defined by the relation 
47ri?^ H pyj = Mconv The quantity tOi is the diffusion velocity of the i— th element at the bottom of the 
convective region that can be expressed as: 



-,5/2 



Ph 



{Ap^i + Vad AT,i 



din P{R^) 
dr 



(33) 



where Ap^i, Aj^i are the diffusion coefficients (see |l4] for details). In the above relation, we neglected, 
for simplicity, the term proportional to the hydrogen concentration gradienlj^ and we took advantage 
of the fact that dlnT/dr = V ad din P/dr at the bottom the convective region, as it is prescribed by 
the Schwartzchild criterion. 



^This term gives a sub-dominant contribution at the present stage and, thus, is neghgible also in the initial phases of 
the evolution when the sun is essentially homogeneous 
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We have to estimate the effect of a generic modification of the SSM on elemental diffusion. In most 
cases, a good description is obtained by assuming that b varies proportionally to the efficiency of 
diffusion in the present sun, i.e. to the the r.h.s of eq.(32) evaluated at the present time. This implies 
that: 

<5A,b = Soji - 5H (34) 

where 5H is the relative variation of the convective envelope effective thickness and Sui is the relative 
variation of the diffusion velocity. By following the calculations described in the appendix B, one is 
able to show that: 



6Dzu = rz6n + rp6P^ 



(35) 



where Ty = 2.05, Tz = 2.73 and Tp = —1.10 and we assumed that all metals have the same diffusion 
velocity as iron. 

By using the above relations and taking into account eqs. ( [25|26| ) one obtains the surface abundances 
variations Alt, and 6Zi, as a function of the free parameters Ai^ni and 5C. Namely, we obtain: 



AYb = Ay AYi^i + Ac5C 
AZb = By AYini + Bc5C 



(36) 



where 



Ay 

Ac 



and 



1 + Dyu 



l + Py.bFiniDybTy 



0.838 
0.033 



(37) 



By 

Be 



Ay 



1-n 

Ac_ 

1-n 



-1.088 
-0.043 



(38) 



Moreover, we can use rels.( 27|2T ) to calculate the variation of the initial metal abundance SZi^i which, 
in our approach, coincides with the variation of metal abundance in the radiative region 5Z(r). We 
obtain: 

6Z{r) = 6Z,^i = Qy AFini + Qc SC (39) 

where: 



Qy = Qo - PY,h Ay {QiTy + Q2Tz) = -0.887 

Qc = -PY,bAc{QlTY + Q2rz)+Tp{Qi + Q2) 



-0.164 



If we had neglected the effect of diffusion, we would have obtained Ay = 0.879 and Ac = 0, By 
and Bc = and Qy = Qo = -1.14 and Qc = 0. 



(40) 
-1.141 
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Figure 5: The coefficients 7/t(?') which define the continuity and the hydrostatic equilibrium equation. 

7 The final set of equations 

We have now all the ingredients to formally define LSMs. The equations obtained in SectJH] namely 
eqs.(30), (36) and (39), relate the chemical composition of the modified sun to the present values of 



the structural parameters 6P{r) and 6T{r), to the energy generation coefficient 5e{r) and to the free 



parameters AYi^i and 5C. These relations can be inserted in eqs.(15) and in the integration conditions 



eqs.( 19 ) and ( 21 ) , obtaining a linear system of ordinary differential equations that completely determines 
the physical and chemical properties of the modified sun. In this section, we give the equations in their 
final form. 

The properties of the sun in the radiative region (r < i?b) are described by: 

[yp 6P + 'jT ST - 6m + -yy Allni + 7e 5e] 
[(7P -1)5P + -fT ST + 5m + A11„i + 7e 5e] 
[p'p 6P + (3'tST-61 + P'y AYini + f3'c6C + 0, 5e] 

[a'p 6P + a'r 6T + 61 + a'y AYi^i + a'c 6C + 6k + a[ 6e\ (41) 
The coefficients jh which define the continuity and the hydrostatic equilibrium equation are given by: 

7P = 1 - Py^p = -1 + PyS_t lY = -Py^Y le = -PY^e (42) 

and are shown in Figj5] as a function of the solar radius. The coefficients and which define the 
transport and the energy generation equation are given by: 



d5m 


1 


dr 




ddP 


1 


dr 


IP 


dSl 


1 


dr 


k 


d5T 


1 


dr 





a' 



Pp + PyCp Pt = Pt + PyCt Py = Priv + PzQy 

ap + dy^p a'rp = Ut + Oty^T Oty = ay£,y + UzQy 



P'c- 



PzQc 



p: 



1 + pY^e 
ae^e (43) 



and are shown in Fig|6] 

The integration conditions at the center of the sun (r = 0) are given by: 



6m 



iPfl SPo + 7T,o + 7yo Al^ni + 7£,o (^eo 
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Figure 6: The coefficients ajj(r) and (3'f^{r) which defines the energy transport and the energy conservation 
equation. 



6P = 
6T = 
SI = 

7t,o = 
-1.142, 



SPo 
STo 



(44) 



-0.254, jY,o = 0.503, 7,,o = 0.161 and P'p^ = 1.647, /3^o = 1.884, 
P'^Q = 0.624. By imposing these conditions, the solution is obtained as a 



where jpfi = 1.167, 
P'co = 0.006, = 

hnear function of the four parameters AYjni, SC, STq, SPq which are fixed by requiring that: 



6m 
6P 
6T 
61 



W-conv 6C 

6C 

A'y AYini + A'c 6C 




at the convective boundary (r = Rh), where Ay = —Py^hAy 
^conv = 0.0192. The solution is, thus, univocally determined. 



0.626, A'(. = -Py^hAc 



(45) 
0.025 and 



The chemical composition of the radiative region can be calculated by using rels.(30) and (39). The 



variation of the density profile 6p{r) in the radiative region can be obtained by the relation: 
6p{r) = -ip{r)6P{r) + -fT{r)6T{r) + 7y(r)A11ni + 7e(r-)(5e(r) . 



(46) 



In the lower layers of the convective envelope the quantities 6P{r), 6T{r), 6p{r) are nearly constant 
(see sectjs]), while the quantity 6m{r) gradually goes to zero, as it is expected being the total mass 
of the sun observationally determined. The surface chemical composition can be determined by using 



rel.(36). We recall that the variations of surface helium and metal abundances are related among each 



other, in such a way that the hydrogen-to-metal surface ratio is unchanged with respect to the SSM. 



Finally, the variation of the convective radius is calculated by using rel.(12). We obtain: 

5i?b = Ty Alini + Tc6C + r« 

where 



(47) 



Ty 



Ay 



-PY,h («^T,b - 4) Ky,t 



1-n 



0.449 
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Ac 

' Cb 
1 



-Py^h (KT,b - 4) + Ky,b 



1-Fb 



«^p,b + 1 
Cb 



-0.117 



-0.085 



(48) 



In the derivation of the above relation, we considered that = 0, 5mb — and we used eqs.(21) and 



(36). 



8 Comparison with full non linear calculations 

In oder to show the vahdity of the proposed approach, we consider four possible modifications of the 
standard input and we compare the results obtained by LSMs with those obtained by the full non linear 
solar model calculations. Three of the considered cases concern with opacity modifications that are 
generally described as: 

K{p,T,X.i)=F{T)K{p,T,Xi) (49) 

where k{p,T, Xi) represent the standard value and F{T) is a suitable function of the temperature. 
Namely, we consider: 

OPAl - Overall rescaling of the opacity by a constant factor F(T) = 1.1. This clearly corresponds to 
introduce the opacity variation: 

SkM = 0.1 



in eqs.(41) for LSM calculations. 



OPA2 - Smooth decrease of the opacity at the center of the sun described by the function: 

F{T) = 1 + ^ 



1 -|- exp 



fTc 



where A = —0.1, Tc = 9.4 x 10^ K, / = 0.1. In our approach this correspond to a variation 6k{r) along 
the solar profile given by: 

6K(r) = 



1 -|- exp 



Tc-T(r) 
fTc 



where T(r) is the temperature profile predicted by the SSM (see left panel of FigjT]). 

OPA3 - Sharp increase of the opacity at the bottom of the convective envelope described by assuming 

F{T) = 1.1 for T < 5 X 10^ K (and F{T) = 1 otherwise). In our approach, this corresponds to: 

5K{r) = 0.1 if T(r)<5xlO*^K 
(5K(r) = otherwise 

The above inequality can be rewritten in terms of a condition on the distance from the center of the 
sun obtaining 5K,{r) = 0.1 for r > OARq (and 6K{r) = otherwise). 



The fourth studied case concerns with modification of energy generation in the sun. More precisely, 
we consider: 
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Figure 7: Left panel: The variation of the opacity profile Sn^r) corresponding to the 0PA2 case. Right Panel: 
The variation energy generation profile Se{r) corresponding to the Spp case. 



Spp - Increase of the astrophysical factor Spp of the p + p d + e'^ + i'e reaction by +10%. In order to 
introduce this effect in LSM calculations, we consider the following variation of the energy generation 
profile: 

Hr) = es,,{r) 6Spp (50) 

where 

/ \ 91ne 

and 5Spp = 0.1. The function 5e{r) obtained in this way is shown in the right panel of FigjTj The 
bump at r ~ O.28i?0 is due to the out-of-equilibrium behaviour of helium-3. 



8.1 Physical and helioseismic properties of the sun 

In Figjsjwe show with solid lines the physical properties of LSMs, obtained by solving the linear system 
of eqs.(|41[), and with dotted lines the results obtained by the "standard" non-linear solar model (SM) 
calculation^^ The four panels corresponds to the input modifications introduced in the previous 
section. We use different colors to show the variations of pressure (red), temperature (black), mass 
(blue) and luminosity (green) as a function of the solar radius. 

We see that a very good agreement exists between LSM and non-linear SM results. To be more 
quantitative, the response of the sun to the input modifications is reproduced at the 10% level or 
better, in all the considered cases. It is important to note that all the relevant features of the 6h{r) 
are obtained, indicating that the major effects are correctly implemented in our approach. The small 
differences between LSMs and non-linear SMs, typically more evident just below the convective region 
and/or at the center of the sun, basically reflects the accuracy of the assumptions which have been 
used to estimate the variation of the chemical composition of the present sun. In the OPAl case, 
a disagreement exists in the 5l{r) behaviour at the center of the sun. This difference, which has no 

^"The term "standard" is used here to refer to the conventional way of calculating solar models, see sectj2] For the sake 
of precision, we use the acronym "SM" to refer to these models. This is done to avoid confusion with the Standard Solar 
Model (SSM) which is intended as our best possible model of the sun, i.e. the model calculated in the conventional way 
and, moreover, by using the best possible choice for the input parameters. 
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Figure 8: Comparison between the physical properties of the sun predicted by LSM (solid lines) and by "standard" 
non-linear SM (dashed lines). In different colors we show the variations of pressure (red), temperature (black), 
mass (blue) and luminosity (green) as a function of the solar radius. 



observable consequences, is mainly due to non linear effects. One should note, in this respect, that the 
constant 10% increase of the radiative opacity induces variations of temperature and pressure which 
are at the ~ 1% level and a much smaller variation of 5l{r) (at the 0.1% level). In this situation, 
cancellations between different first-order competing contributions may occur. 

In TabjT] we present the variations of the initial abundances and of the photospheric chemical 



composition obtained by using rels.(39) and (36). Moreover, we show the variation of the convective 



radius 5R\j calculated according to eq.(47). Finally, in FigjD we present the variations of the density 



profile 6p{r) calculated according to eq.(46) and the variation of the squared isothermal sound speed 
u = Pj p given by: 

5u{r) = 6P{r) - 5p{r) . (52) 

We see that LSMs reproduces the results obtained by "standard" calculations with a very good accuracy. 
In particular, the variations of the helioseismic observables are obtained within 10% or better (unless 
they are extremely small). All this shows that our approach is sufficiently accurate to use LSMs as tool 
to investigate the origin of the present discrepancy with helioseismic data. In a separated paper, we 
will use LSMs to analyze the role of opacity and metals in the sun [26] . 
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Figure 9: Comparison between the variations of density (red) and of squared isothermal sound speed (black) 
predicted by LSM (solid lines) and by "standard" non-linear SM (dashed lines). 



8.2 Neutrino Fluxes 

As a final application, we calculate the solar neutrino fluxes predicted by LSMs and we compare them 
with those obtained by using "standard" non-linear SM calculations. The flux of solar neutrinos 
can be expressed as: 

^y = jy2 j dr r'^ p{r) ny{r) (53) 

where D is the sun-earth distance, ny{r) is the total number of neutrinos produced per unit time and 
unit mass in the sun and the index v = pp, B, Be, N, O labels the neutrino producing reaction according 
to commonly adopted notations. If we expand to first order, the relative variation of the fiux can 
be expressed as: 

5(^y = ^ J dr {r/D)'^p{r)n^{r) X 

[n^,p{r) 5p{r) + ny;T{r) 5T{r)n„^Y{r) AY{r) + ni,,z{r) 5Z{r) + n^,spp{r) 6Spp] (54) 



where is the SSM value and: 



din 



dlnp 



+ 1 



SSM 
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OPAl 


0PA2 


0PA3 


Spp 




SM 


LSM 


SM 


LSM 


SM 


LSM 


SM 


LSM 




0.016 


0.017 


-0.0056 


-0.0058 


0.0021 


0.0019 


0.0015 


0.0016 




-0.018 


-0.016 


0.000 


-0.001 


-0.012 


-0.012 


-0.010 


-0.011 




0.014 


0.014 


-0.0037 


-0.0036 


0.0038 


0.0038 


0.0031 


0.0034 


5Z^ 


-0.018 


-0.018 


0.0049 


0.0047 


-0.0050 


-0.0049 


-0.0040 


-0.0044 


SRb 


-0.0020 


-0.0020 


-0.0067 


-0.0070 


-0.014 


-0.015 


-0.0058 


-0.0064 




-0.011 


-0.010 


0.0045 


0.0052 


-0.0020 


-0.0011 


0.0090 


0.0092 




0.13 


0.13 


-0.067 


-0.064 


0.017 


0.016 


-0.11 


-0.11 




0.27 


0.27 


-0.17 


-0.17 


0.029 


0.028 


-0.27 


-0.28 




0.14 


0.14 


-0.10 


-0.094 


0.003 


0.004 


-0.21 


-0.22 


S^o 


0.21 


0.22 


-0.14 


-0.14 


0.012 


0.012 


-0.29 


-0.31 



Table 1: Comparison between the predictions of LSM and "standard" non-linear SM for the initial and surface 
chemical abundances, the convective radius and the solar neutrino fluxes. Note that the absolute variations are 
reported for Helium, whereas the relative variations are shown for all the other quantities. 



ny,z{r) 



din rii^ 



dlnT 

din riu 



dlnZ 
d In n,. 



dY 
din rin 



SSM 



SSM 



SSM 



dlnS, 



pp 



SSM 



(55) 



The derivatives nj/j- (r) have been calculated numerically by assuming that the abundances of secondary 
elements (helium-3, carbon-12 and nitrogen-14) can be estimated as it is described in appendix A. The 
symbol |ssM indicates that we calculate the derivatives along the density, temperature and chemical 
composition profiles predicted by the SSM. 

The term nj/^Sppl?") ^•S'pp is introduced to describe the effects of a variation of S'pp on the rates of 
neutrino producing reactions. It clearly holds npp^SppC'^) = 1) since the pp— neutrinos are produced by 
p + p ^ (i + e^ + i/g. Boron and Beryllium neutrinos arc influenced by S'pp, because this parameter 
determines the hclium-3 production rate (through deuterium) and, thus, also the rate of ^He-|-^He — 
Be -|- 7 reaction. As a consequence, we have: 



"'Be,Spp(r) = nB,Spp( 



dlnS, 



pp 



(56) 



SSM 



where is the helium-3 abundance. The CN- cycle efficiency, instead, does not depend on the value 
of S'pp and, thus, we have nN,Spp(?') = ''^o,pp('') = 0. 
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Figure 10: The functions 4>u.j{r) defined by eq.(58) 



It is useful to recast eq.(54) in the form: 
S'^u -- 



dr [(A.,p(r) 5p{r) + (P.^rir) 5T{r) + (l,,y{r) Ay(r) + (A.,z(r) 6Z{r) + (/).,spp(r) 5S, 



ppJ 



where: 



J dr p(r) ny{r) 



(57) 
(58) 



The functions (pujif) are displayed in Fig, 10 They show explicitly the well-known fact that the boron, 
beryllium and CNO neutrinos are produced in a more internal region with respect to pp-neutrinos and 
that they more strongly depend on temperature. The complicated behaviour of the functions ^Nj{r) 
at r ~ O.15i20 is due to the out-of-equilibrium behaviour of carbon-12 abundance. The A'^— neutrinos 
originate, in fact, from the decay of ^'^N, which is produced by ^^C + p N + 7 reaction. The 
sun was born with a relatively large amount of carbon-12, which has been converted by CN-cycle into 
nitrogen-14 in the more internal regions. The ^^C abundance is, thus, larger where CN-cycle is less 
effective and, as a consequence, a non negligible production of A^— neutrinos occurs at relatively large 
radii, where equilibrium conditions do not hold. 

In Tab{l] we compare the LSM results, i.e. the values 6^,y obtained by applying eq.(p7|) to the 
solutions of eqs.(41), with the results obtained from "standard" calculations. We see that an excellent 
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agreement exists, for all the cases considered in this paper All this show that our approach can 
be used as a tool to investigate the dependence of the solar neutrino fluxes on the input parameters 
adopted in SSM construction. 

9 Conclusions 

We have proposed a new approach to study the properties of the sun which is based on the following 
points: 

i ) We have considered small variations of the physical and chemical properties of the sun with respect 
to Standard Solar Model predictions and we have linearized the stellar equilibrium equations to relate 
them to the properties of the solar plasma (see sects. 
a ) We have derived simple relations which allow to estimate the (variation of) the present solar compo- 
sition from the (variation of) the nuclear reaction rates and elemental diffusion efficiency in the present 
sun (see sectjG]). 

As a final result, we have obtained a linear system of ordinary differential equations (see sect j7| which 
can be easily solved and that completely determine the physical and chemical properties of the "per- 
turbed" sun. 

In order to show the validity of our approach, we have considered four possible modifications of 
the input parameters (opacity and energy generation profiles) and we have compared the results of our 
Linear Solar Models with those obtained by "standard" methods for solar model calculations (see sectjs]). 
A very good agreement is achieved for all the structural parameters (mass, luminosity, temperature, 
pressure, etc.) and for all the helioseismic and solar neutrino observables. 

We believe that our approach can complement the traditional method for solar model calculations, 
allowing to investigate in a more efficient and transparent way the role of the different parameters and 
assumptions. In particular, it could be useful to study the origin of the present discrepancy between 
SSM results and helioseismic data. 
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^^It is remarkable the fact that all the neutrino fluxes are correctly reproduced, even in the cases in which some 
compensations occur. See e.g. the OPA3 case, in which a very small value for 5$]v flux is obtained (with respect to the 
other fluxes), as a consequence of the peculiar behaviour of the ^N,j{T) functions. 
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A: The abundance of secondary elements 



Secondary elements are those elements which are both created and destroyed in a reaction chain or in a 
reaction cycle. The relevant secondary elements to understand energy and neutrino production in the 
sun are ^He, ^^C and ^^N. 



Al: The '^He abundance 



The evolution of helium-3 abundance, X3, is described by the equatior 



dt 



-{av) 



pp 



9 



(av) 



33 



X3Y 
12 



(av) 



34 



(59) 



where {crv) pp is the reaction rate per particle pair of the p + p ^ d + e'^ + Ve reaction, while {(tv)^^, and 
((7^)34 refer to ^He +^ He -^'^ He + 2p and ^He +^ He -^'^ Be + 7 reactions, respectively. In the most 
internal region (r < O.25-R0), the rates of these reactions are fast with respect to the sun evolutionary 
times and equilibrium is achieved: 



3,oq 



3r M34 
8 (cjt;)33 



-1 + 




In the outer regions {R > 0.2Rq), we can neglect the contribution from ^He +^ He 
obtain: 

^ ^ ^2 



dt 



C3 - V3 Xi 



(60) 

Be + 7 and we 
(61) 



where C3 = (3 p (erf )pp)/(2 m„) and P3 



(p(o-i')33)/(3m„) 

1 - 



If we assume that the factors P3 e 



C3 are approximately constant over the time scale T3 = C3 ^ during which ^He is produced, the above 



equation can be explicitly solved. We obtain: 



X3 — X- 



3,eq 



tanh 



^3. 



eq 



(62) 



where X3,max = C3 1 represents the helium-3 abundance which would have been produced in the inte- 
gration time t if we had neglected the ^He destruction processes (and we implicitly assumed that the 
helium-3 initial abundance is negligible). 

In this work, we need to calculate the helium-3 response to a generic modification of the solar 
properties. We use the following approach. We calculate the ^He equilibrium value in a generic solar 



model by using eq.(60). We assume that eq.(|62j) is valid at each point of the sun and we estimate the 
value X3^max in the SSM by inverting it (and by using the SSM-abundance X3). Then, we calculate 
the value of X3_niax in the modified solar model by considering that X3 max scales as -'^3,max oc p Sp 



^pp- 



Finally, we calculate the value of X3 in the modified sun by using eq.(62) with the modified values for 
X3^eq and X3^max- The obtained results reproduce with very good accuracy the variation of the ^He 
abundances in all the cases considered in this paper. 



^We assume that deuterium is at equilibrium, i.e. deuterium production and destruction rates are equal. 
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A2: The ^^C and ^^N abundances 

The evolution of carbon- 12 and nitrogen- 14 abundances, indicated with X12 and X14 respectively, is 
described by the equations^ 



dX 



12 



dt 

dXu 

dt 



12 p 

ruu 
Up 

rUu 



XX 



14 



14 
X X12 



12 



(0-u)l,14 



XX 



12 



12 

X Xii 



14 



(0-f)l,12 
(0-l')l,14 



(63) 



where (c7v)i,i2 and (o"?j)i^i4 are the reaction rates per particle pair of "^"^C+p — N-I-7 and ""^^N-l-p -^^^ 
O -I-7, respectively. It is useful to rewrite these equations in terms of the variables ry = X14/I4 — X12/I2 
and N = X12/I2 + X14/I4. We obtain: 



drj 

di 
dN 

~dt 



Cr^N-Vr,r] 





(64) 



where = p X{{av)i^i2 — (o"u)i,i4)/m'u and 2?^ = p X{{av)i^i4, + (cf )i,i2)/w-u- From the above 
equations, we see that is constant and that the equilibrium value for rj is given by: 

(o-y)i,i2 - (o-t^)i,i4 



(o-y)i,i2 + {crv) 



1,14 



(65) 



In the assumption that the coefficients and T>ri are approximately constant over the time scale 
Tr^ = l/{CrjN), the solution for r] can be explicitly calculated, obtaining 



V = Veq + {VO - Veq) exp 



^eq 



(66) 



where 770 is the initial value and r/max = C N t represents the r/ value which would have been obtained 
in the integration time t if we had neglected the "destruction" term in eq. ( 64 ) . 



In analogy to what was done for helium-3, we use the following approach to calculate the response 
of carbon-12 and nitrogen-14 to a generic modification of the solar properties. We calculate the 7?eq 
value in a generic solar model by using eq.(65). We assume that eq.(66) is valid at each point of 
the sun and we estimate the value r/j^^x ™ the SSM by inverting it (and using the SSM-values fj 
and r/g). Then, we calculate the value of r/max in the modified sun by considering that it scales as 
^max oc N pX {{crv)i^i2 — (o'v)i,i4). Finally, we calculate rj in the modified solar model by using eq.(66) 
with the new values for r/eq, ??max and 770- The carbon-12 and nitrogen-14 abundances can be calculated 
from r] and A'^ by the simple relations: 



X12 = 6iN-rj) 
Xi4 = 7{N + rj) 



(67) 



The obtained results reproduce with very good accuracy the results obtained by full numerical calcu- 
lations in all the cases considered in this paper. 



^We neglect the NO-cycle, which is largely sub-dominant with respect to the CN-cycle in the sun. 
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B: Derivation of eq.(35) 



The effect of elemental diffusion on the convective envelope chemical composition can be estimated 
from: 

(5A,b = SuJi - 6H (68) 

where 6H is the relative variation of the convective envelope effective thickness, while 6uJi is the relative 
variation of the diffusion velocity of the i— element at the bottom of the convective region. One should 
note that the quantity 5L0i is defined as: 



(69) 



and, thus, involves the difference between the diffusion velocities evaluated at two different points. 
By taking into account the expression for cjj given in eq.(33) and by considering that din P/dr ~ 
-{G^Mq/RI)(p/P) at r = flb, we obtain: 



5uji 



5 r(i?b) - r(i?b) 

2 



P(i?b) - P(i?b) 



where A 



tot,i 



r(i?b) P(i?b) 

Ap^i + Vad^T,i- This can be rewritten as: 
5 



25i?b + ^^^Alb 
aY 



(70) 



5uji 



-5rb-5Pb-25i?b + ^^^^Ayb 



5 91nr(i?b) 
2 



5 In ; 



91nP(Pb) 
dlnr 



6R^ 



(71) 



where, according to the notation convention adopted in this paper, we have indicated with 5Tb = 
[T(Ry,) -T(Rk)]/T(Rk) and (5Pb = [P(Rh) -P(Rh)]/P(Rh)S We note that the last term in the above 
equation can be neglected since f?lnT/f?lnP = Vad — 2/5 at the bottom of the convective region. 

The effective thickness of the convective envelope is defined by the condition H = Mconv/(4vri?^/>b). 
The relative variation 6H is thus given by 



6H = 6Mc. 



9AR X 91np(Pb) 

2dRb - Oph ORb 

omr 



(72) 



If we consider that Mconv 
region, we obtain: 



dr 47rr^/9(r) and that 5p{r) is approximately constant in the convective 



6Mc, 



Sph - = SRh 
li 



By using this relation and taking into account eqs.(72) and (71), we arrive at the expression: 

5 „„ dlnA 



5D 



6n - + 



totji 



dY 



Ayb + 



(9 In p(i?b 
dlnr 



Rh 

H 



6Rh 



(73) 



(74) 



Finally, by using the expression for 5Riy given in eq.(12) (with dmy^ = and (5/b = 0) and by taking 
advantage of the integration conditions at the bottom of the convective region expressed by eqs.(21), 
we obtain the final relations: 



SDY^h = Ty STi, + Tp6Pk + Tf,6K^ 
5Dzh = Tz 5T\, + Tp5Ph + Ti^6Hih 



(75) 



*We assumed that dlnT{R^,)/dlnr ~ dlnT(Ri,) /dlnr and dlnP{Rb)/dlnr ~ dlnP{Rb)/dlnr. 
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where Ty = 2.05, Tz = 2.73, Tp = —1.10 and F^; = —0.06 and we assumed that all metals have the 
same diffusion velocity as iron. In the calculations presented in this paper, we neglect for simplicity 
the terms proportional to 6k\j which generally give a very small contribution. 
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